Conserved Spin and Orbital Angular Momentum Hall Current in a Two-Dimensional 
Electron System with Rashba and Dresselhaus Spin-orbit Coupling 
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We study theoretically the spin and orbital angular momentum (OAM) Hall effect in a high 
mobility two-dimensional electron system with Rashba and Dresselhuas spin-orbit coupling by in- 
troducing both the spin and OAM torque corrections, respectively, to the spin and OAM currents. 
We find that when both bands are occupied, the spin Hall conductivity is still a constant (i.e., 
independent of the carrier density) which, however, has an opposite sign to the previous value. The 
spin Hall conductivity in general would not be cancelled by the OAM Hall conductivity. The OAM 
Hall conductivity is also independent of the carrier density but depends on the strength ratio of 
the Rashba to Dresselhaus spin-orbit coupling, suggesting that one can manipulate the total Hall 
current through tuning the Rashba coupling by a gate voltage. We note that in a pure Rashba 
system, though the spin Hall conductivity is exactly cancelled by the OAM Hall conductivity due to 
the angular momentum conservation, the spin Hall effect could still manifest itself as nonzero mag- 
netization Hall current and finite magnetization at the sample edges because the magnetic dipole 
moment associated with the spin of an electron is twice as large as that of the OAM. We also eval- 
uate the electric field-induced OAM and discuss the origin of the OAM Hall current. Finally, we 
find that the spin and OAM Hall conductivities are closely related to the Berry vector (or gauge) 
potential. 

PACS numbers: 73.43.Cd, 73.63.Hs, 75.47.-m, 85.75.-d 
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I. INTRODUCTION 



Spin transport electronics or spintronics in semicon- 
ductors has become a very active research field in con- 
densed matter mainly because of its potential applica- 
tions in information storage and processing and other 
electronic technologies pj and also because of many fun- 
damental questions on the physics of electron spin 0- 
Spin current generation is an important issue in the 
emerging spintronics. Recent proposals of the intrinsic 
spin Hall effect are therefore remarkable 0, Q. In the 
spin Hall effect, a transverse spin current is generated in 
response to an electric field in a metal with relativistic 
electron interaction (spin-orbit coupling). This effect has 
been considered to arise extrinsically, i.e., by impurity 
scattering [j| . The scattering becomes spin-dependent in 
the presence of spin-orbit coupling, and this gives rise to 
the spin Hall effect. In the recent proposals, in contrast, 
the spin Hall effect can arise intrinsically in hole-doped 
(p-type) bulk semiconductors Q and also in electron- 
doped (n-type) semiconductor heterostructures Q| due to 
intrinsic spin-orbit coupling in the band structure. This 
intrinsic spin Hall effect offers an exciting possibility of 
pure electric driven spintronics in semiconductors, where 
spin-orbit coupling is relatively strong and which can be 
more readily integrated with well-developed semiconduc- 
tor electronics. 

A large number of theoretical papers have been writ- 
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ten addressing various issues about the intrinsic spin Hall 
effect. In 0, a systematic semi-classical theory of spin 
transport is presented, resolving a discrepancy between 
the prediction of Q and the Kubo formula result. In 
0, an orbital-angular- momentum (OAM) Hall current 
is predicted to exist in response to an electric field and is 
found to cancel exactly the spin Hall current in the spin 
Hall effect. In ||, however, ab inito relativistic band 
structure calculations show that the OAM Hall conduc- 
tivity in hole-doped semiconductors is one order of mag- 
nitude smaller than the spin Hall conductivity, indicating 
no cancellation between the spin and OAM Hall effects 
in bulk semiconductors because of the orbital quench- 
ing by the cubic crystalline anisotropy. There is also an 
intensive debate about whether the intrinsic spin Hall ef- 
fect remains valid beyond the ballistic transport regime 
0, firi . 1 1 1 1 ] . On the other hand, experimental measure- 
ments of large spin Hall effects for the Rashba two di- 
mensional electron gas and for n-type bulk semiconduc- 
tors have just been reported |12lll3j ]. although more work 
is needed to firmly establish the intrinsic or extrinsic na- 
ture of the results. 

At present, an urgent current issue in spintronics re- 
search is about the a ppr opriate definition of the spin cur- 
rent 0, 0, 0, ITU Il8j. In almost all previous studies 
of the intrinsic spin Hall effect, the spin current is intu- 
itively defined as the expectation value of the spin and 
velocity operators, namely, ^{v, s z }, where {, } is the an- 
ticommutator defined by {A, B} = AB + BA. Similarly, 
the OAM current is defined as the expectation value of 
^{v, L z }. Here, s z and L z are the z-component of spin 
and OAM operators, respectively, and v is the veloc- 
ity operator. However, this conventionally defined spin 
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(OAM) current is not conserved in systems with spin- 
orbit interaction |(J. Consequently, many fundamental 
questions on the intrinsic spin Hall effect in semicon- 
ductors remain unresolved. Very recently, in Reflla a 
proper definition of the conserved effective spin current 
is established for systems with spin-orbit coupling, and 
the conserved spin current density is defined as 



(1) 



where J s = Re\W'i{v, s z }^] is conventional spin cur- 
rent density and P r is the torque dipole density which 
arise from spin torque T" = Re[^' ^[s z , Hq]^] where 
Hq is the Hamiltonian of the system. The commutator 
[, ] is defined by [A, B] = AB - BA. The effective spin 
current density J s then satisfies the standard continu- 
ity equation [16j of ^f- + V • J a = where the spin 
density is defined by S z = ty^s^. A derivation of the 
effective spin continuity equation in Rashba-Dresselhaus 
system is given in Appendix A. Within this definition of 
the effective spin current, the unphysical intrinsic spin 
Hall effect in insulators with localized orbitals vanishes, 
and the Onsager relation between the spin Hall effect 
and inverse spin Hall effect is ensured. [l6( Furthermore, 
this new definition of the spin current predicts opposite 
signs of spin Hall coefficients for a couple of semicon- 
ductor models such as the Rashba and fc-cubed Rashba 
Hamiltonians. |l6( 

In this paper, we extend the theory proposed in Ref. 
ITfilto two-dimensional electron systems with both Rashba 
and Dresselhaus spin-orbit coupling. Furthermore, we in- 
troduce the concept of the effective OAM current by in- 
cluding an OAM torque correction term, and investigate 
the OAM Hall effect in both pure Rashba system and sys- 
tems with Rashba-Dresselhaus spin-orbit coupling. We 
derive the effective OAM continuity equation for the 
Rashba-Dresselhaus system (Appendix B). Also in this 
paper, we argue that in a pure Rashba system, though 
the OAM Hall conductivity is found to exactly cancel 
the spin Hall conductivity, there nevertheless would be 
nonzero magnetization current and finite magnetization 
at the sample edges because the magnetic dipole moment 
associated with spin angular moment is twice as large as 
that of the OAM. Finally, we find that there are interest- 
ing relations between the Berry vector potential and the 
spin and OAM Hall conductivities. 



II. MODEL HAMILTONIAN AND LINEAR 
RESPONSE CALCULATION 



which stems from the structural inversion asymmetry. 
Here, o = (a x ,a y ) and a z are the three Pauli matrices. 
The other is the Dresselhaus SO coupling which results 
from the bulk inversion asymmetry, if the heterostructure 
is made of semiconductors without spatial inversion sym- 
metry such as semiconductors in the zincblende struc- 
ture [2£| ■ The Dresselhaus term is given by 



(3) 



Therefore, the full Hamiltonian for the 2DEG with 
Rashba-Dresselhaus SO coupling can be written as 



H = 



P 

2m 



H R + H L 



(4) 



where m is the effective mass of the electrons in the 
2DEG. Interestingly, though the Dresselhaus coupling co- 
efficient [3 is fixed for a given structure, the Rashba cou- 
pling str eng th A can be tuned by a gate voltage by up 
to 50 % |2l| . thereby providing an opportunity to study 
the interesting interplay between both types of the SO 
coupling. The Hamiltonian has been solved exactly by 
several authors. The eigenstates can be written as 



where n = ±1 is band index and 



-i0 (k) 



tan0(k) 



Xky f3k x 



(5) 



(6) 



(7) 



The corresponding eigenenergies is given by 



Zm 

where k = (k Xl k y ), k = |k|, tan(0) = k y /k x , cos6>(k) = 
(Acos0— (3 sin (f>)/ j ((f)), sin#(k) = (Asin0— j3 cos 4>)/^ ((f)) 
and j(<f>) = y/(X 2 + (3 2 ) - 2A/3 sin(2</>). It can be proved 
that the difference of two Fermi momenta hk F is given 
by (kp —kp) = 2m'y((j))/h 2 . The band structure consists 
of two energy bands which are degenerate at the centre 
of the 2D momentum space k = 0, as schematically illus- 
trated in the left panel of Fig. 1(a). The spins associated 
with the eigenstates all lie in the xy plane, as shown in 
the right panel of Fig. 1(a) and in the left panel of Fig. 
1(b). 



A. Rashba-Dresselhaus Hamiltonian 

For a two-dimensional electron gas (2DEG) confined in 
a semiconductor heterostructure, two major spin-orbit 
(SO) interaction terms are usually present. One is the 
Rashba term 0] , 



(P x e z ), 



(2) 



B. Spin Hall conductivity 

Let us consider an uniform electric field E applied in 
the y-direction. The total Hamiltonian in this case is 
given by H = Hq + eEy, where — e is the electron charge. 
Let us treat the term eEy as a small perturbation. To 
determine the spin transport coefficient, we start with 
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FIG. 1: (Color online) (a) The schematic bandstructure of the 
Rashba-Dresselhaus spin-orbit coupling system with |A| 7^ \f3\. 
The outer sheet is for n = +1 and the inner sheet, for n = — 1. 
The two sheets touch at k = 0. If A = f3, the sheets be- 
come degenerate at line k x — k y (see, also, Fig. 1 in [23^. 
The right panel is the top view of the two bands. The ar- 
rows labeled on the two bands denote the spin directions of 
the associated fc-points, respectively, (b) Initially, the spins 
are in the xy-p\a,ne (the left panel). At a short time to af- 
ter an external electric field E (indicated as the large hor- 
izontal (red) arrow) is applied along the y axis, the bands 
move in the -k v direction with the distance &k y = eEto/h 
(see Subsec. II. E). Each spinor feels an effective magnetic 
field and will do precession. For k x > 0, the spins tend 
to tilt down and for k x < 0, the spins tend to tilt up 
(the right panel). The z component of the spin is given by 
$(sz) = -eEh(X 2 - (3 2 )cos(f>/4nk 2 -y 3 ((f)) (c) Before an exter- 
nal electric field is applied, all the eigenstates carry zero or- 
bital angular momentum (OAM), i.e., {L z )o = (see Subsec. 
II.E.). When the electric field is turned on, the z component 
of the OAM is induced as shown in the right panel (small 
vertical arrows) . A linear response calculation would give the 
z component as: 8{L Z ) — eEfi(X 2 — (3 2 ) 2 cos cf> / 4nk 2 j 5 ((f)) . 



the linear response Kubo formula in the clean limit. The 
conventional spin Hall conductivity is [22| 



where j s x z = l/2{v x , s z } is the usual spin current op- 
erator, v(k) = dH (k)/hdk, H (k) = e - lkx iJ e lkx , 
s z = fcr z , |nk) is the eigenstate of the nth band with 
momentum k, and / n k is the Fermi-Dirac distribution. 
The index n ^ n' denotes no intraband transition. Sim- 
ilarly, the spin torque response coefficient can also be 
obtained Hfl: 



X(q) = -7T Yl^ nk ~ /«'k+q] 

n^n' k 

(nk\r s (k, q)|n'k + q) (n'k + q| v(k, q)|nk) 

(E„\t — -En'k+q) 2 



(9) 



where r fl (k,q) = |[r s (k) + r s (k + q)], v(k,q) = ±[v(k) + 
v(k + q)] with r s = j^[s z ,H ). In the above equation, 
a s x ° y is the transport coefficient of conventional spin cur- 
rent and a^y can be determined from the dc response 
of the spin torque dipole: P r = i?e[iV q x(q)]| q ^o • E, 
namely, crjj = Re[idxy(q)/dq x ]\ q ^o [l6|. Here x(q) is 
the response coefficient of spin torque to an electric field 
with finite wavevector q. Before doing summation in 
the above Kubo formula, we expand the term inside the 
summation to the first order of q. At zero temperature, 
/n'k+q ^ /n'k + q • (dE n , k /dk)df n > k /dE n / k , where the 
first term is a step function while the second term was 
involved with a Dirac delta function. We assume that 
the two bands are all occupied by the electrons and the 
Fermi energy is larger than the SO splitting. The spin 
Hall conductivity is 



'xy: 



(10) 



which is the linear response of the conserved effective 
spin current J s to the electric field ^(|. It is straight- 
forward to calculate the spin Hall conductivity (the de- 
tailed derivation is given in Appendix C). We find that 
a x°y — —sign(X 2 — P 2 )e/8ir is the conventional spin Hall 
conductivity and a r x a y = sign(X 2 — f3 2 )e/4ir is the con- 
ductivity due to the spin torque correction. Therefore, 
the spin Hall conductivity defined as the response of the 
conserved effective spin current is: 



\ 2 =I3 2 



0. 



(ii) 



which is independent on the carrier density. Remarkably, 
with the spin torque correction term a xy included, the 
sign of the effective spin Hall conductivity is opposite to 
the conventional spin Hall conductivity reported in |23j 
and H. 



y ^ t ^ ' [/nk /n'k] 
n^n' k 

Im[(nk\j s x *\n'k) (n'k|« v |nk)] 



K'k) 2 



(8) 



C. Orbital angular momentum Hall conductivity 

The usual OAM current operator was introduced as 
1 



:{v x ,L z }, 



(12) 
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where L z is the z component of the OAM operator. In 
the Rashba Hamiltonian and with this definition of the 
OAM current, the spin Hall effect is always accompa- 
nied by the OAM Hall effect which was first noted by 
Zhang and Yang [7J. In a SO system, such a definition 
has the same problems as the conventional spin current 
operator because the OAM current is not a conserved 
quantity. Therefore, we derive the continuity equation 
for the effective OAM current (see Appendix B) . We find 
^§f- + V • J Q = T 2 °, where C z = ^L z ^> is the z compo- 
nent of the OAM density, J Q = Re[^'f°H] is the OAM 
current density and j 00 = ±{v, L z }. T° = Re[^r ^} 
is the torque density, where t — i [L z , H ] . In ana- 
logue to the conserved effective spin current operator 
proposed in ^(|, the average torque in the bulk is zero, 
and hence we have y J dVT° = 0. The torque den- 
sity can be written as the divergence of torque dipole 
density P°(x), i.e., T z = -V ■ P?(x). Substituting 
T° = — V • P° into continuity equation of OAM, we ob- 
tain J a = 3 a + P° as the effective OAM current. The 
torque dipole density vanishes outside the bulk, and we 
can write f v dVP° = f v dV '(-xV • P?) = J v dVxT z °. 
Therefore, the unique form of torque dipole density is 
P° = Re^' (xt )^]. In short, we may define an effective 
OAM current operator: 



Jo 



~dt ' 



(13) 



It has an extra term xr D which is the correction term due 
to the OAM torque. The corresponding OAM current 
density J Q — i?e^ / ^(x) l j 5'(x) defined as the expectation 
value of this current operator satisfies the standard con- 
tinuity equation of ^if- + V • J Q = 0. As for the spin 
Hall conductivity jig , the OAM Hall conductivity has 
two parts: 



'xyl 



(14) 



where the first term is the usual OAM Hall conductivity 
and the second term comes from the OAM torque correc- 
tion. Thus, the OAM Hall conductivity can be calculated 
from Eq. JSJ) and Eq. J^J) by substituting j** and t s with 
j°° and t = i [L z , H ] , respectively. It can be shown 
(Appendix C) that the OAM Hall conductivity defined 
as the linear response of the effective OAM current is: 



A 2 +/T e \2 ^ o2 
^2g¥' A >P 

0, A 2 - f3 2 



A 2 - 



which is a function of the ratio A//3 and is also indepen- 
dent of the carrier density. 



ket of system will travel on a close path and return to 
initial state ket after a period. The final state ket must 
coincide with initial state vector, apart from a phase fac- 
tor. Berry [25| has shown that the state ket will acquire 
an additional phase factor as the system undergoes the 
adiabatic evolution. The phase factor accompanying the 
adiabatic evolution is called Berry phase. Berry phase 
cannot be removed by any gauge transformation for a 
closed path [26J. In solids, the Bloch state also acquire 
a Berry phase if the applied perturbation can make a 
constraint such that k travels adiabatically on a closed 
path in the Brillouin zone j2?J • The Berry phase of Bloch 
states \rik) for a closed path C can be expressed by 



A n (k)-dl, 



(16) 



where A„(k) = (nk|(— i)^|nk) is defined as Berry vec- 
tor potential (or Berry connection). Berry curvature is 
defined as 0„ = Vk x A n . The Berry vector poten- 
tial and curvature are the salient characteristic of energy 
band structure and hence have important applications 
in transport properties of carriers. It has been shown 
that the equation of motion of Bloch electron has an ex- 
tra anomalous velocity in terms of Berry curvature of 
Bloch states |2^. In bulk p-type semiconductors with 
spin-orbit coupling, the existence of fc-space magnetic 
monopole (Berry curvature) in the degeneracy T point of 
band structure results in a transverse force exerting on 
spin Q . The connection between the dissipationless spin 
Hall conductivity and Berry curvature has been shown 
in [29j. Here we further show that in two-dimensional 
Rashba-Dresselhaus systems, we can write the spin Hall 
and OAM Hall conductivities in terms of Berry curvature 
and vector potential. 

With the eigenstates |nk) given by Eq. JjjJ, we can 
prove that the matrix element (n'k|i;j,|nk) satisfies the 
following relation (for n' ^ n): 

(n'k\v y \nk) = ~(£ nk - £ n , k )(n'k|H) — \nk) (17) 

where (n'k|(— i)-^-\nk) = —\-§^~- The Berry vector po- 
tential is then given by 



A=(nk| H )|| nk) = _l|. 



(18) 



(15) Substituting Eas. itTTl . and (O into Eq. @ and noting 



that 



t s « = -cr So 

xy yx 



cf,~ I30L we have 



eh 



n k 



/nk 



(19) 



D. Berry vector potential and Hall conductivity 



where w„k is defined as w„k = {E ri \^— E-n^) / h. Similarly, 
we also obtain, 



When the Hamiltonian of a physical system is param- 
eterized by periodically changing environment, the state 



n k 



'xy 



mV 
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where a°° y = -a°° x was used The spin Hall and OAM 
Hall conductivities are therefore related to Berry vector 
potential. The Berry curvature is fl(k) z = [Vk x A] z = 

— |^~^| nS^ (k). The corresponding Berry phase is 
then given by the path integral of the Berry vector poten- 

A 2 /3 2 

tial, i.e., $ = § c d\- A= J R Q,(k) z dk x dk y = - \ X A_pi\ n = 
—sign(\ 2 — /3 2 )7r, where i? is the region containing the 
origin in the k-plane. The Berry curvature £l z is zero for 
k 7^ but non-zero for k = where the two bands of 
the Rashba-Dresselhaus system are degenerate. As for 
the spin transverse force for the Luttinger Hamiltonian 
given in Ref. l29i the corresponding spin transverse force 
in Rashba-Dresselhaus system can also be given in terms 
of the Berry curvature. Even though the Berry curvature 
vanishes for k ^ 0, the spin Hall current can nonetheless 
occur due to the Aharonov-Bohm-like effect at k = [3lJ . 

Let us now discuss the effects of the different choice of 
the eigenstates on the above findings. In, e.g., Ref. [32|. 
the eigenstates (|nk)) 

= 71 ( "T ) ^ = 71 U-0 • < 21 > 

were used for the Rashba-Dresselhaus system. We find 
that the unitary transformation defined by the matrix 

U--( ~ l + el9 (22) 

with UU^ = Wll = 1, will transform the eigenstates of 
Eq. © to that of Eq. J2TJ), i.e., \nk) = U\nk). Using 
the eigenstates in Eq. (|2*T)> . the Berry vector potential is 
given by 

,~ , , . s d , ~. 1 89 

A + --( + k|H)-| + k)^--- 

A -" ( - k|( - l) 9k | - k)=+ 2ak 

that depends on the band index n. We may define A„ = 
nA for both the bases of JSJ and i|21|l . As a result, the 
Berry phase given by Eq. 1)16(1 also depends on the band 
index n [^3, i.e., 

$± = j> A± (k) • dl = Tsign(\ 2 - (3 2 )n = ±$, (24) 

where $ is the Berry phase when the basis of Eq. (JSJ is 
used. In contrast, the signs of the Berry phase $ for the 
two bands are the same in the choice of the eigenstates 
of Eq. 10. Nevertheless, we can write $„ = ?i<i> with 
n = ±1. Note that the charge Hall effect would be zero 
when the basis of Eq. I|21|) is used, and would be nonzero 
if the basis of Eq. © is chosen, as have been pointed out 
in Ref. IH 

Using A = — 1# (Eq. IT5)l . we can show that A = 
[k x A],a , where a = ft) = % and [k x A] 2 = 

— (A 2 — (3 2 )/2j 2 (4>) which depends only on (j>. After some 



vector algebra calculations, the spin Hall conductivity in 
Eq. (|19|l can be written as 

< = ^E£<"-»[>'xA] I ao ¥ f^, (25) 

where (hp — kp) = ?H^W was usec |_ Because {k F — 
kp n ) = n(k^ — kp), Eq. becomes 

<% = E / dl ■ n t k x A "]- a ° ( 26 ) 

where A„ = nA with n = ±1. Using the definition of 
Berry phase in Eq. I)16|l . the spin Hall conductivity can 
be written as 

<v = 77^2 E f dl ■ (" A «) ^E'*- ( 2? ) 

n JC n 

where $„ = n<i> with n = ±1, in agreement with Ref. 
l23l In other words, the spin Hall conductivity can be 
written as a line integral of the Berry vector potential or 
a surface integral of the Berry curvature. Note that Eq. 
(12711 holds irrespective of the choice of the eigenstates. 

For the OAM Hall conductivity, similarly, Eq. l(2Tj|l can 
be rewritten as 

a °° y = ^E£ dl -[ kxA «]* a °- ( 28 ) 

This expression has also been given in Ref. l33l where 
[Hk x A] 2 = -H(X 2 - /3 2 )/2 r y 2 ((j>) was regarded as the 
OAM of the eigenstates in the absence of any applied 
electric field. We believe that this interpretation is in- 
appropriate because the expectation value of the conven- 
tional OAM operator depends on the choice of bases, as 
will be shown in next Subsec. below. In next Subsec, 
we also argue that before an electric field is applied, the 
OAM of the eigenstates is zero. Finally, we also find 
that the ratio of the spin to OAM Hall conductivity is 

-Iff = 2[k x A] z |0 =m7r where m is an integer. 



E. Origin of the orbital angular momentum Hall 
current 

All electrons have an intrinsic spin of ^. Before an 
in-plane electric field is applied, the spins of the elec- 
trons are all aligned in the icy-plane, as shown in Fig. 
l(a)-(b). When an in-plane electric field is applied, the 
SO coupling gives rise to not only a spin transverse force 
on a moving electron 0, |3f| but also an effective SO 
magnetic torque. The Rashba-Dresselhaus Hamiltonian 
can be written as (Hji(k) + Hjj(k)) — a ■ B e // where 
~B e ff = A(k x e z ) — P(k x i x — k y e y ) is the effective SO 
magnetic field. The dynamics of the z component of spin 
can be derived from the Heisenberg equation of motion 
and we obtain da z (t)/dt = Jf -[a z ,^ + H D + H R + eEy] = 
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±-[a z ,H D + H R ] = ^[a z ,ai]Bif f . Using the commu- 
tation relation of spin matrix [<7j,<7 7 -] = lit^y-u^, we 
have (ds z (t)/dt) = -|(s(t) x B e //) z , where s = §<?. 
^From the equation of motion of electron in fc-space, 
we also have Sk y = (k y (to) — k y (t = 0)) = — an d 
Sk x — 0. Therefore, after a short time to, the Fermi 
surface (i.e., the circle for n = ±1) would move along 
-k y direction with the distance eEto/H. This implies 
that the variation of the effective SO magnetic field is 
STi e ff = (Xe x + (3e y )dk y . Consequently, each spin feels 
the effective SO magnetic torque s(0) x ^B e // and tilts 
out of the xy plane. The quantum dynamical analysis of 
spin for the Rashba-Dresselhaus system has been given 
in p^ |. Here we use the quantum perturbation method 
to evaluate the response quantities. 

Let us expand the wave function to first order of elec- 
tric field, |nk)' = Ink^ + lnk) 1 ^) where the perturbed wave 
function is 



(29) 



The expectation value of the z component of spin can 
be evaluated by '(nk|s 2 |nk)' = (s z )o + S(s z ) where 
(s z )o = (nk|s z |nk) = and the z component of the 
spin in first order of electric field is given by S(s z ) = 
2i?e(nk|s z |nk)W = ~eEh(X 2 - (3 2 ) cos <?!>/4nfc 2 7 3 (</>)• 
Therefore, for k x > 0, the spins on the outer (inner) 
sheet tend to tilt down (up) and for k x < 0, the spins on 
the outer (inner) sheet tend to tilt up (down), as shown 
in the right panel of Fig. 1(b). This results in trans- 
verse spin Hall currents with spin polarization in the z 
direction. 

Now consider the OAM L z = H(x xk) z . The ex- 
pectation value ' (n\i\L z \rik)' = (L z )o + S(L Z ) where 
(L z ) = (nk\L z \nk) and 5{L Z ) = 2i?e(?ik| J L z |?ik)( 1 ) . In 
contrast to the spin case, before an electric field is ap- 
plied, as argued below, all the eigenstates carry zero or- 
bital angular momentum, i.e., (L z )q = 0, as illustrated 
in the left panel of Fig. 1 (c) . In the absence of the ap- 
plied electric field, the diagonal matrix elements of the 
OAM depend on the choice of the eigenstates, though the 
off-diagonal matrix elements of the OAM do not. For 
example, if we choose the eigenstates of Eq. (5), to- 
gether with the conventional position operator x = ?Vk, 
(L z ) = (nk|x x hk\nk) = -h{X 2 - ^]/2 7 2 (^), a vame 

which has also been obtained in Ref. Ha. On the other 

/ e -tf(k)/2' 

hand, if we use the eigenstates |nk) = I 

(L z )o — 0. This is of course unsatisfactory, and therefore, 
we propose to define a gauge invariant position operator 
X = iVk -I- A(k) to resolve the problem (see Appendix 
D). We can show that with this gauge invariant position 
operator X, (L z )o — 0, irrespective of the choice of the 
phase factor of the eigenstates. Therefore, we believe 
that all the eigenstates of the Rashba-Dresselhaus sys- 
tem carry zero OAM in the absence of applied electric 
fields. Importantly, the other quantities such as S(s z ) 



and 8(L Z ), which contain only the n ^ n' matrix ele- 
ments, are independent of the choice of the phase factor 
and also the choice of the position operator. 

As an in-plane electric field is applied, the dynamics 
of the z-component of the OAM can be obtained by 
means of Heisenberg equation of motion: dL z (t)/dt = 
|(/ix x k) z = -|(s(t) x B' eff ) z + (x x F) z where 
= Be//|A->_A and F = -eE = -eEe y . The second 

term (x x F) z = ftx x k is the classical external torque 
which causes the orbital motion of the electrons and de- 
pends on the choice of the origin of the coordinate system. 
The first term — -|(s x B^) z is the effective SO magnetic 
torque. If we use the gauge invariant position operator, 
the classical torque (x x F) z does not contribute to the 
variation of the OAM in a short time tg. We can show 
that |(SXxk) z =-§(s(i) xB' eff ) z -X{t)eE + k y n z , 

where D, z — Vk x A ~ <5' 2 ' (k) and is zero because 
k 7^ 0. After a short time to, (X(0))Sk y vanishes be- 
cause (X(0)) = and the z component of the OAM for 
each eigenstate is induced by the effective SO magnetic 
torque and external torque only. In our calculations, 
we treat the potential eEy as a perturbation and use 
again the quantum perturbation method instead of solv- 
ing the Heisenberg equation of motion. We find that 
the z component of the OAM is S(L Z ) = eEh(\ 2 — 
( 9 2 ) 2 cos^/4nfc 2 7 5 W = S(s z )[{\ 2 - P 2 )h 2 {<j>)}, as illus- 
trated in the right panel of Fig. 1(c). Therefore, under 
the in-plane electric field along the y axis, the electrons 
on the outer (inner) sheet with k x > drift toward the 
x axis, carrying the downward (upward) tilted spinors 
as well as finite positive (negative) OAM S(L Z ), and the 
electrons on the outer (inner) sheet with k x < drift to- 
ward the x axis, carrying the upward (downward) tilted 
spinors as well as finite negative (positive) OAM S(L Z ). 
This gives rise to the OAM Hall current. It is inter- 
esting to note that in the pure Rashba system ((3 = 
and (A 2 - (l 2 )h 2 (<t>) = 1), the S(s z ) and S(L Z ) have the 
same magnitude but the opposite signs and hence cancel 
each other exactly. This is due to the fact that the z- 
component of the total angular momentum is conserved 
in pure Rashba system. 



III. SYSTEMS WITH BOTH RASHBA AND 
DRESSELHAUS SPIN-ORBIT COUPLINGS 

The calculated total and decomposed spin and OAM 
Hall conductivities are summarized in Table I. The to- 
tal spin and OAM Hall conductivities are displayed as a 
function of the ratio \X/f3\ in Fig. 2. The total angular 
momentum Hall conductivity (er*° + a^ s y + a° x ° y + a x ° y ) is 
(see Table I) 



xv ~\ 0, A 2 ^/3 2 



(30) 



It is clear that in 2DEG systems with the Rashba- 
Dresselhaus SO coupling, the total angular momentum 



7 



TABLE I: Total and decomposed spin (a) and orbital angu- 
lar momentum (b) Hall conductivities of a two-dimensional 
electron system with Rashba-Dresselhuas spin-orbit coupling. 

(a) Spin Hall conductivity 





A 2 > P 2 


A 2 = f3 2 


A 2 < f3 2 


a s ° 


e 





e 


u xy 


s» 




T s 


e 





e 


xy 


4-7T 


4tt 


s 


e 





e 


® xy 


Sir 


8tt 



(b) OAM Hall conductivity 





A 2 > f3 2 


A 2 = /3 2 


A 2 < 1 


u xy 


\' A +P 2 e 





\ I +0'- t e 


\\' 2 -P 2 \ Srr 




T Q 

® xy 


e 





\ £ +P z e 


\X i —0'\ 4tt 


A 2 -,3 2 | 4ir 


a 

^ xy 


V-h9* e 





V+P* e 




A^-^l 8tt 



Hall conductivity is in general not zero, except that 
A 2 = 1 or (3 — (pure Rashba SO coupling). The 
sign of the total angular momentum Hall conductivity is 
always negative, and, as will be discussed shortly, this is 
because the Hall conductivity is dominated by the nega- 
tive OAM Hall conductivity. Furthermore, when the two 
SO coupling strengths are comparable, the total angu- 
lar momentum Hall conductivity is very large, suggest- 
ing the interesting possibility of tuning the angular mo- 
mentum Hall effect by varying the Rashba SO coupling 
strength. This large total Hall conductivity results from 
the large OAM Hall conductivity in both the conven- 
tional and present definitions of spin current (see Table 
I). The conventional OAM Hall conductivity is the same 
as that given in Ref. |3iJ This singularly large OAM con- 
ductivity in the region that |A//3| approaches to unity, 
could lead to spontaneous magnetization, as suggested 
in Ref. |3j Nevertheless, we believe that this singular 
behavior of the OAM conductivity near A ~ (3 is unphysi- 
cal and is perhaps due to the neglect of the disorder in the 
present calculations. The disorder due to, e.g., impurity 
scattering and electron-electron interaction is known to 
have pronounced effects on the spin Hall effect [sL IToLfTTj . 
The infinite large OAM conductivity is expected to be 
suppressed by the disorder effects in real 2DEG systems. 
Nonetheless, further calculations taking into the disorder 
effects are beyond the scope of the present paper. 

Shen recently pointed out that in the 2DEG sys- 
tems with both Rashba and Dresselhaus couplings, the 
spin current along the z direction is antisymmetric with 
respect to an unitary transformation: <j x — > a y ; <j y — > a x ; 
o z — > — a z . This antisymmetry makes the conventional 
spin Hall conductivity changes sign at A 2 = (3 2 . It is 
interesting to note that the spin torque Hall conduc- 
tivity and hence the total spin Hall conductivity also 
change sign at A 2 = (3 2 when one moves from the re- 
gion where Dresselhaus coupling dominates to the region 
where Rashba coupling dominates (Table I and Fig. 2). 
This shows that our calculated spin torque and total spin 
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FIG. 2: (Color online) Spin {a% y , solid line) and orbital angu- 
lar momentum {a% y , dashed line) Hall conductivities versus 
the strength ratio |A//3| of the Rashba (A) to Dresselhaus (/?) 
spin-orbit coupling. 



Hall conductivities also obey the requirement of this an- 
tisymmetry. We also note that as for the case of pure 
Rashba coupling ^||, the spin torque Hall conductivity 
is twice as large as the conventional spin Hall conductiv- 
ity but has an opposite sign (Table I), giving rise to the 
result that the conserved spin Hall conductivity has the 
same size but opposite sign to the conventional spin Hall 
conductivity. 

On the other hand, all the (total, torque and conven- 
tional) OAM Hall conductivities do not change sign at 
A 2 = (3 2 , as shown in Table I and Fig. 2. As for the case 
of spin Hall effect, the torque OAM Hall conductivity is 
twice as large as the conventional OAM Hall conductivity 
but has an opposite sign, resulting in that the effective 
OAM Hall conductivity has the same size but opposite 
sign to the conventional OAM Hall conductivity (Table 
I). The conventional OAM Hall conductivity has been re- 
ported in Ref. .33, and our result is consistent with this 
previous calculation. 



IV. SYSTEMS WITH EITHER PURE RASHBA 
TERM OR PURE DRESSELHAUS TERM 

A. Pure Rashba spin-orbit coupling 

In the case of j3 = and A ^ 0, the Hamiltonian re- 
duces to the pure Rashba Hamiltonian. The spin torque 
operator becomes 



T ? = h Sz P+H R ] 

in Lrn 



and the y-component of velocity operator is 



v y = ~k ax 

u m n 



(31) 



(32) 
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The torque correction term for the spin Hall conductiv- 
ity is e/47r, and for the OAM Hall conductivity is — e/47r, 
which exactly cancel each other. This result is confirmed 
by the realtion [L z , Hr] = —[s z , Hr] that is present in 
the Rashba Hamiltonian. Including the torque correction 
term, we get an effective spin Hall conductivity of e/87r, 
and an effective OAM Hall conductivity of — e/87r. This 
implies that when the torque correction term is taken into 
account, the spin Hall current is still exactly canceled by 
the accompanied OAM Hall current and there is no total 
angular momentum current in the pure Rashba system. 
This interesting observation has been reported before in 
Ref. in the context of the conventional definition of the 
spin and OAM currents. This may be expected because 
[s z + L Z , Hr] = and the conservation of the total angu- 
lar momentum therefore must be obeyed no matter what 
definitions of the effective angular momentum currents 
one adopts. 

Nonetheless, we want to point out here that the spin 
Hall effect in the pure Rashba Hamiltonian can still man- 
ifest itself and be detected in several ways, even though 
the total angular momentum Hall current is zero. It is 
well known that the magnetic dipole moment associated 
with the spin of an electron is fj, a = — 2s/is and the one 
associated with the OAM of an electron is fj, = — L^b 
(see, e.g., [3^). Consequently, although the total angular 
momentum (s z + L z ) current is zero, the total magneti- 
zation [(jitg + fig) = — {2s z + L z )\xb\ current would not 
be zero and will give rise to finite magnetization at the 
edge of the sample. Therefore, the spin Hall effect in the 
pure Rashba Hamiltonian can in principle be probed in at 
least two ways. As noted in Ref. 0, one is by measuring 
the electric field induced by the nonzero magnetization 
current |3?1 ] and the other is the magnetization at the 
sample edges. Recently, it is suggested in Ref. O that 
spin current can be detected by measuring time depen- 
dent magnetization precession. 



B. Pure Dresselhaus spin-orbit coupling 

In the case of A = and (3^0, the Hamiltonian 
reduces to the pure Dresselhaus Hamiltonian. In this 
case, the total, torque correction and conventional spin 
and OAM Hall conductivities have the same signs (Ta- 
ble I). For example, the total spin Hall conductivity is 
— e/87r, being the same as the total OAM Hall conduc- 
tivity, thereby giving rise to a total angular momentum 
Hall conductivity of — e/4ir. In the Dresselhaus Hamilto- 
nian, the total angular momentum is not conserved, i.e., 
[s z + L z ,Hb>] = —2i/3(a x p y + a v p x ) ^ 0. Instead, we find 
[s z , Hd] — [L z , He] or [s z — L z , Ho] — 0. Therefore, the 
spin and OAM Hall conductivities would add up rather 
than cancel each other, in contrast to the Rashba Hamil- 
tonian. 

We notice that the spin Hall conductivity in the clean 
limit is constant in either pure Rashba system, or pure 
Dresselhaus system, or the mixture Rashba-Dresselhaus 



system [see Table 1(a)]. Interestingly, in Refs. ITol and 
UlT it is pointed that the constant spin Hall conductiv- 
ity in pure Rashba system would result in an unphys- 
ically growth of the magnetization with time. This is 
due to the fact that the conventionally defined spin Hall 
current in the pure Rashba system turns out to be pro- 
portional to the time derivative of the magnetizaion of 
the system. This linear relation between the conven- 
tional spin Hall current and the time derivative of the 
magnetization should also exist in the pure Dresselhaus 
system because the two systems are related by the uni- 
tary transformation: a x — > a y and a y — ► a x , as men- 
tioned above in Sec. III. Indeed, we find that the con- 
ventional spin current in the pure Dresselhaus system is 
3x" = ji v x, cr z } = ^cr z p Xj where the relation {a t , aj} = 
2Sij was used. Furthermore, it can be shown that the 
spin precession and conventional spin current satisfy the 
re l ations: f^M = §{p x a z )(t) = ^(f x *(t)) and 

d<<T * t (t)> = ^-{jy'{t)). Therefore, the constant conven- 
tional spin current in the clean limit in the pure Dres- 
selhaus system would also lead to the unphysical con- 
sequence that both (a y (t)) and (<J x (t)) depend linearly 
on time, resulting in an infinite growth of the magneti- 
zation. In contrast, because the conserved spin current 
satisfies the continuity equation [Eq. (1)], we find no such 
relation that the conserved spin current is proportional 
to the time derivative of the spin operator, thereby, free 
from the artifact discussed above, as also pointed out in 
Ref. HI 

V. CONCLUSIONS 

In conclusion, we have calculated the spin Hall coduc- 
tivity in the Rashba-Dresselhaus Hamiltonian with the 
spin torque correction in the absence of disorder, and 
find it to have the same magnitude but an opposite sign 
to the result reported before 0, |2^| . The spin Hall con- 
ductivity in the absence of disorder is still a constant in 
pure Dresselhaus ( or Rashba ) system even when the 
torque correction is considered. We also introduce the 
conserved effective OAM current and find that in gen- 
eral, the OAM Hall effect does not cancel the spin Hall 
effect in the 2DEG with the Rashba-Dresselhaus spin- 
orbit coupling. The OAM Hall conductivity depends sig- 
nificantly on the strength ratio of the Rashba to Dressel- 
haus spin-orbit coupling, suggesting that one can manip- 
ulate the total Hall current through tuning the Rashba 
coupling by a gate voltage [3^, l4fj | . We argue that in a 
pure Rashba system, though the spin Hall conductivity 
is exactly cancelled by the OAM Hall conductivity due to 
the angular momentum conservation, the spin Hall effect 
still manifest itself as nonzero magnetization Hall current 
and finite magnetization at the sample edges because the 
magnetic dipole moment associated with the spin of an 
electron is twice as large as that of the OAM. We show 
that the spin and OAM Hall conductivities have a simple 
relation to the Berry vector (or gauge) potential. We also 
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calculate the electric field-induced OAM and discuss the 
origin of the OAM Hall current. 
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APPENDIX A: THE CONTINUITY EQUATION 
FOR THE EFFECTIVE SPIN CURRENT 



In this appendix, we derive the continuity equation 
for the effective spin current. A similar derivation can 
also be found in |15|. The definition of spin density is 
S z = ^s z ^, where s z = \a z . Using ih-^^S = and 
H = Hq + eE • x, wc obtain the following equation: 

at in 



p 2 

i -*) t ^*] (Al) 



1 

ih'~ i '2m v 2m 
+ [&a z H R y - (H R ^a z V] 

+ [&a z H D V - (H D ^a z ^]}, 

where we have used [s z ,x] = 0. First, we combine the 
following two equations: 

(* T (Tz^) = ( ) cr z \P + a z ^ 

dxi dxi dxi dxi dxi dxi 



d d& a _ ( _d__d_yu a y . d g 

dxi dxi z dxi dxi z dxi dxi z 

Using the definition of momentum operator p = —ihV 
and (* t pcr 2 \I') t = -(p<jt)(7 2 Yl/ 5 W e obtain : 

p • [2Re{¥pa z ^)j = [¥p 2 a z ^ - (p 2 * t )<r z *]. (A2) 

Next, consider the term [^a z H R ^ - {H R ^a z % which 
contains only the Rashba Hamiltonian. We obtain the 
following equations: 



[W<t,Hr9 - (H R yya z y} 
n n 



and 



[¥a z H R V - {H R ^a z ^\ 

= ^[* f P ' ^ - jrP ■ [^a z (a x e z )*] f . 



(A4) 



We combine Eq. (|A3|) and Eq. (|A4|) : 

[&*,H R 9 - (H R ^a z y} 

2i\ A 
= — i?e[* f (p • - -p • Re[Wa z (a x e 2 )*]- 

Now consider the term : [^<j z H D ^- (H D ^a z ^] which 
contains only the Dresselhaus Hamiltonian. From the 
second term (Hp^y <J Z ^ ', we obtain: 

{H D 0a z ^ = - ~[-p x &<T x + p y ^a y ]a z ^ 


= ~Px{Wct x (t z -$) - ^&p x a x a z V 
n n 

- ~p y (&a y a z y) + jr^p y (jya z y 
= ~ j:{p x {^<JzV x y) ~ Py(^a z a y ^)} 



-^\p x a x - Py(Jy)<7 z ^. 



(A6) 



and then: 

-^\p x a y +p y cr x )^ 



-Py[*Wi/*]} (A7) 

-2i0 



-Re[¥(p x (Ty+pya x )^} 



+ j:{p x Re[¥a z a x ^] - p y Re[¥ a z a y ^}}, 

where we have used the relation: [9'a z Hn^ — 
(H D yy<T z S>] = -[&o- x H D V - (H D yy <r x Vy* . Substi- 
tuting Eqs. (|X2|) . ljA5j) and (|AT|) into IjAljl . we obtain 



= - V • Re[^a z {— - -a X e z - ^(a x e x - a y e y ))^\ 
m n 



O \ OR 

^i?e[* t (p • - -^Re[tf(p x a y + p y a x )^]. 



(A8) 



Consider the commutator [a z , Hq], and we have 



1 r P 2 rr rr 1 2A - / \ 

(A9) 

The velocity operator in the Rashba-Drcsselhaus system 
satisfy the following relation: 



v = l[x,i/] = l[x,ff ] 
in in 

P A /3 
= -{a x e z j - —{a x e x - <J y e y ). 

TT\t lb lb 



(A10) 
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Finally, substituting Eq. l|A10f> and Eq. <|A9|) into Eq. 
H A8[l . we obtain the continuity equation of the effective 
spin current 

— [¥s z m] = -V • Re[¥'f^] + Re[^T s ^}, (All) 

where we have used s z = |ct z , = ^{v, s z } and 
t s = jjr[s z ,Ho]. With the definition of spin density 
S z — 9's z *, spin current J s = i?e[ v E' 1 'j SzV E'] and spin 
torque 7 Z = Re^^T^], Eq. (|A11|) can be rewritten as 



dt 



V • J s = T/ 



(A12) 



which is in agreement with Ref. [if: 



APPENDIX B: THE CONTINUITY EQUATION 
FOR THE EFFECTIVE OAM CURRENT 

In this appendix, we derive the continuity equation 
for the effective OAM current. Let us consider the z- 
component of the OAM density C z = ^L z * and its 
partial time-derivative ih-^C z . Applying ih-j^* = H* 
where H = Hq + eE • x, we have 

= ^L Z H ^- (H Q ^L z y + ¥[L z ,e~E-x.}y 

2 2 



^[L^eE-x]* 



(Bl) 



First, using the commutation relations [p, L z ] — ihe z x p 
and [L Zj p 2 ] = 0, we can show that: 

{¥L z p 2 y-(p 2 &)L z *} 
= p-[4' t pL z *-(p* t )(L z *)] 

= p [¥{p,L z }^} -p 2 {¥L z ^)+p- [¥iM z x p*]. 



(B2) 



For the Rashba Hamiltonian H R , we have 

(H R ^L Z * = ^[-Py^a x +p x ^a y }L z ^ 
= ^p- [&(a x e z )L z *] 
+ ^ t (p y (T x -p x a y )L z y 



and 



{¥l z H r V - (H R ^L Z *} 

A + (B3) 
= &[L» H R ]V - -p • [¥(a x e z )L z *}. 



For the Dresselhaus Hamiltonian Hd, we obtain 

{&L,H D V - (H D ^L Z *} 

B (B4) 
- &[L Z ,H D ]V - ■ [&(a x e x - a y e y )L z V}. 

Substituting Eqs. JEU, CEU and (JHH into Eq. JED, w e 
have 

ih^- t (¥L z v) 

2m 

- r^O? x e z )L z * - ~&(<T x e x - a y e y )L z y] (B5) 

+ &[L Z , H R + H D ]* - -*-p 2 (*%#) 
2m 

+ y~ V [^ihe z x p*]. 
Putting Eq. I|A10|) into Eq. l|B5)l . we have 

= - V • [*ti{ v , L z }#] + *t I eE • x]* 

+ ^^[L z , Ho}* - ^V 2 (*+L,*) - p • [tftg, x ptf] 

Finally, it can be shown that p • [^e z x p*] is purely 
imaginary, and ^L z * — Re(^L z *) is real for both 
the eigenstates of the Rashba-Dresselhauls Hamiltonian 
and the Bloch states ^k. The C z continuity equation 
becomes : 



dt 



(# t £ z 1') = -V-i2e[* t j OQ *]+i2e[* t r *]+-Re[* t T B *] 

(B6) 

where j°° = \{v, L z } , r G = ^[Z z ,if ] and t e = 
-^[L z ,eEi ■ x]. With the OAM current density J D = 
Re[^j°°*], the OAM torque density 7° = i?e[#V #] 
and the classical torque density T Z E = fle[#V B #], Eq. 
(IB6I) can be rewritten as 



dC z 
dt 



V-J = T Z ° + T Z 1 



(B7) 



which is the OAM continuity equation. 

The torque 7 E — Re[WTE*] has a classical ana- 
logue. It can be regarded as the rotational torque mo- 
ment due to the force eE exerted on a particle located 
at the position x with respect to the origin of the coor- 
dinate system. This can be seen from the commutator 
[L z , eE ■ x] = —ihe(x x E) z , and 



T* = -ite[(#Tx#) x eE] 



(B8) 



^From the space integral / dV7 z E = —Re[(J dV^x*) x 
eE] z , we obtain / dV7 z E = -i?e[(#|x|#) x (eE)] z , where 
/ dV^x* = (*|x|*). Expanding the (*|x|*) in pow- 
ers of electric field, we have (#1x1*) = (# | x l#o) + 
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(x)W + o(E 2 ), where |^o) satisfies the unperturbed wave order of q x : 
equation ihMA^Qj = Hq^q). Therefore, we obtain 

JdVT z E = -Re[(V \x\* Q )x(eV)} z +o(E 2 ). However, we (nk|r(k, q)|n k + q)(n k + q| W (k, q)>k) 

should demand that the expectation value of the OAM in _ _ n ^ j : ®® 

the unperturbed system is zero. This would imply that m y dk x 



the expectation value of the position operator must be A/3 

zero in the absence of the external electric field, as dis- h~ y + ^ y cos ^^ + ^ x sin ( 2 ^) 

cussed in Appendix D. Finally, we obtain that the OAM qq „ 

continuity equation + <lx-^-{K cos(26») - k y sin(20)) + sin(20)] 

OK x Z / C3 x 

2 

- ! -r[k x + k v sin(26>) - k x cos(26») 

^ + V • 3 = T° (B9) go 

dt +q x — {k y cos(26) + k x sm(28)) 

which is valid in first order of the electric field. "^2 cos(20))] + o(q x ). 



APPENDIX C: SPIN AND OAM TORQUE HALL 
CONDUCTIVITY 



In this appendix, we derive the spin and orbital angu- 
lar momentum (OAM) torque Hall conductivities (<r xy ) 
in the Rashba-Dresselhaus system. The conductivity in 
pure Rashba system can be obtained by setting = in 
the A 2 > 2 condition. The total torque conductivity is 
given by 



We also need to expand the (/ nk - fn'k+q) / (E n k 
-En'k+q) 2 to first order of q x . We have 



/nk /n'k+q /nk /n'k 



(Enk — Sn'k+q) 2 (E n k — E n > k ) 2 

9i?n'k ro /nk — /n'k 



q x - 



dk x (E nk - E n , k f 
o(q x ). 



(E nk - E n , k f 



(C4) 



Assuming the Fermi energy is lager than the spin-orbit 
splitting and using using (kp~ — hp) — 2m^((f>)/h 2 and 

a xy = — Im[Xy^], we can write the torque conductivity 



vly = {° T x V + <?l° y ) = Re[ id Xy(fl)/dqx\ qx ^o], (CI) 



where Xyio) is defined as 



Xy(q) = ^ X] ^2ifnk - /n'k+q] 
n^n' k 

(nk|r(k, q)|n'k + q) (n'k + q\v y (k, q)|nk) 



(Enk — i?n'k+q) z 



(C2) 



be rewritten as: o T xy 



Eq. (IHTl) . 



To evaluate Xy m Eq. (|C2|) , we first expand the matrix 
elements (nk|r(k, q)| n'k+q) (n'k+q|u(k, q) y |nk) to first 



2tv 



G(A, /?,</>) 
7(0) 



,vr, (C5) 

e/3 f 271 ,,H(X, 0, <j>) cos tf> 



Btt 2 Jo 



where G(A, 0, </>) and H(X, 0, <f>) are given by 

(/? 2 -A 2 )sin 2 0(Asin(20)-/3) 



G(A, /?,</>) 



(7(0)) s 



H(\,0,4>) = A sin (f> + cos <j> 



where r(k,q) = |[r(k) + r(k + q)], v(k,q) = ±[v(k) + 
v(k + q)] and r = j^[s z + L Z , H ] is the total torque. The 
equation (jClfl suggests that we can choose q = (q x , 0, 0) 
for simplicity. By the choice of q = (q x ,0,0), the X B (q) 

can be expanded in power of q x , i.e., x a (q) = XiP + 
Xy^Qx + Xy 2 \ x + °{l x )i an( i the torque conductivity will 

Re[ixy ] = — Im[xy^] by using i(<t>) 2 7(<W 



2 2 (A sin - cos 0) cos(20) 
+ (A ~ P) » 
/ A 2 | ^ 2 ^( Acos '? !) + / ?sin '? !) ) sin(20) 



(7(^)) 2 



-2A/3 



(Acos0 + fismi 
(7(^)) 2 



Using [ 7 (0)] 2 = (A 2 + /? 2 ) - 2A/3sin(2</>), we find that 



H cos < 



can be written as: 



r G(A, /3, ^) 1 ff(A, 0, 4>\ , = (/3 2 _ A 2 } _2A sin cos - 



7(0) 2 7 (^)2 



7(0)4 



(C6) 
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^From the above formula, we can easily check that Eq. 
(|C6[) the is zero when we set A 2 = (3 2 . From Eq. 
(|C6j) . it follows that the torque conductivity vanishes in 
pure Rashba system (i.e. (3 — 0) because the integral 
J d(f> sin (j> cos </> is zero. Our next step is to work out the 
integrals by using Residue method [4l]]. The crucial in- 
tegrals are the followings (|A| > \f3\): 



d4> 



2tt(A 2 +/3 2 ) 
(A 2 - (3 2 ) 3 
• _ 2ir\f3 
(A 2 - (3 2 f 



(C7) 



APPENDIX D: GAUGE INVARIANT POSITION 
OPERATOR 

In this appendix, we shall consider a method to ob- 
tain the gauge invariant position operator and resolve 
the problem that the expectation value of the OAM de- 
pends on the choice of the eigenstates in the absence of 
applied electric field. In crystalline environment , th e po- 



sition operator cannot be simply set to x = i-^ |42J. We 
find that the OAM will depend on the choice of the eigen- 

For example, 



states when we use the definition: x 



"die 



the following states can also be chosen as the eigenstates 



For |/3 1 > |A|, we can exchange the Rashba and Dressel- 
haus couplings A < — > (3 in the integrals. Substituting 
Eq. |C6| and Eq. ifOfjl into Eq. (Hf>) . we obtain the 
torque conductivity: 

- ' ^ (C8) 



2tt |A 2 - (3 2 



Since the commutation relation [s z + L z , + Hr] = 0, 
the Rashba coupling disappears in numerator of a xy . For 
r s fl ee ±-[s z ,H R (k)] = X((J x k x + (Jyky), we can use the 
same steps as the calculation of <J xy . The new H S (X, (3, <p) 
and G s (\,(3,(j>) can be obtained by changing the A and 
/3 in H(X,/3,<j>) and G{X,(3,<f>). We have 



2tt 



+ 



eX 



16tt 2 



7(0) 

2 \ 1M H S {X,(3,4>) cos0 

2 ^ G(A^A^) 
7W 



(C9) 



H{X <-> (3,4>) cos0 

' hW 2 



A 2 



Att\X 2 -P 2 \' 
Therefore, the spin torque Hall conductivity is 

e f3 2 e A 2 . 



■ xy 



(~ 4tt \X 2 -(3 2 \ ' 4tt \X 2 -f3 2 \ J 
e X 2 - (3 2 . 2 2 e 



(CIO) 



and the OAM torque Hall conductivity is 

f3 2 e X 2 



'xy 



(- 



4tt|A 2 -/? 2 | 4tt|A 2 -/3 2 | 



) 



e X 2 +f3 2 
~4tt \X 2 -/3 2 \' 



(Cll) 



For the pure Rashba system, we can take (3 = in the 
A 2 > (3 2 condition. For the pure Dresselhaus system, we 
can take A = in the A 2 < (3 2 condition. 



|nk; M) 



; iM(k) 

~7T~ 



me 



i8(k) 



(Dl) 



where M(k) is the phase factor and n = ±1. It can be 
proved that the eigenstates of Eq. (|D1|) satisfy the re- 
lation Hq = J2 n Enk\nk; M)(nk; M\ where the eigenen- 
ergy is defined in Eq. J7J. Using the eigenstates 
of Eq. (|D1(1 . the matrix element (rik; M\L z \rik; M) is 
zero if M = —9/2 but nonzero if M = —8. Fur- 
thermore, if the eigenstates of Eq. I|21|l are chosen, 
the signs of (+k|L z |+k) and (— k|L 2 |— k) will be dif- 
ferent. Let us now consider the phase transformation 
of |nk) = e lA "( k )|nk) where the phase A„(k) depends 
on band index n. We note that the eigenstates of Eq. 
(J5J will be transformed to that of Eq. I|21|) by the ap- 
propriate choice of the phase A„. Actually, we have 
|+k) = (-i)| + k) and |-k) = e w \ - k). Let us also 
introduce the operator X' ee x + A' where x = i-J^ and 
A' is the vector potential. Using the above gauge trans- 
formation [2(1, we obtain 



(nk|X>k) = (rak|x|nk) + A' - 
Consider the specific transformation: 

A' = / 9A « 



OAn 

9k ' 



dk 



(D2) 



(D3) 



Eq. HD2() becomes the gauge invariant form 43] : 

(nk|X'|nk) = (nk|X|nk), (D4) 

where the operator X is defined as X = i-^ + A. We find 
that A' satisfies the transformation (|D3(I if A' is defined 
as 



A; = (nk|(-i)^|nk) 



(D5) 



which is the Berry vector potential when the eigenstates 
|nk) are used. The vector potential A in the eigenstates 
|nk) can be defined as A„ = (nk|(— i)J^\nk) which gen- 
erally depends on band index n. Therefore, the gauge 
invariant position operator X will depend on band index 
n: 



X ± =4 + A ± . 



(D6) 
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We find that the expectation value (nk|X„|nk) van- 
ishes: (nk|X„|?ik) = 0. It can also be shown that 
(n' 7^ nk|X„|nk) = {p! ^ nk|i^|nk) by the use of 
(n'k|nk) = S n / n where 5 n i n is equal to unity for n = n' 
and zero for n =/= n' . In terms of the gauge invariant po- 
sition operator, the OAM operator (x x k) 2 is replaced 
by (X„ x k) z . As a result, the expectation value of the 
OAM operator (L)o = h(nk\X n \nk) x k is zero for all the 
eigenstates in the absence of the applied electric field. In- 
terestingly, note that the OAM operator does not always 
depend on band index n. For example, if the eigenstates 
of Eq. © are used, the Berry vector potential is inde- 
pendent of n, and so is the gauge invariant operator X„. 

Let us now show that the calculated OAM Hall con- 
ductivities remain unchanged no matter whether the 
conventional x or gauge invariant X„ position opera- 
tor is used. First of all, the gauge invariant velocity 
operator ^-[X-t, Ho(k)] is the same as the original one 
^-[x, Ho(k)] because the vector potential A± commutes 
with -ffo(k), namely, [A±, i/ (k)] = 0. The gauge in- 
variant OAM torque operator ^-[(X± x k) z ,Ho] is also 



the same as the conventional one ^j[(x x k) z ,iJo] be- 
cause (A± x k) z commutes with H (k). Nevertheless, 
the OAM current operator ^{v x ,L z } has an extra term 
of (A± x k) z v x when the gauge invariant position op- 
erator X± is used. For the eigenstates of Eq. (|D1(I . 
we find that (nk;M\v x \ - nk; M)(-nk; M\v y \nk; M) = 
-[(A 2 + /3 2 ) sin 9 cos 6 + \(3}/h 2 and A± = |(2^ + 
gg) are all purely real [44), On the other hand, Eq. 
® needs only the imaginary part of the matrix ele- 
ments (A± x k) z (nk; M\v x \ -nk; M) (-nk; M\v y \nk; M). 
This shows that the extra term (A± x k) z v x does not 
contribute to the Kubo formula. In general, using 
|nk) = e lA ™|nk), we have (rtklt^ln'k) (7i'k|wj,|nk) = 
(nk\v x \n'k) (n'k\v y \nk) . With the eigenstates of Eq. 
(I21[) . we obtain that (nk\v x \—nk)(—nk\v y \nk) = — [(A 2 + 
/3 2 )sin6>cos6» + \(3]/h 2 and A± = T^f^, which are all 
purely real. Therefore, we conclude that the OAM Hall 
conductivities a xy and a xy remain the same even if we 
replace x with X+ or X_ . 
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